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We reveal that the periodic radius modulation peculiar to one-dimensional (ID) peanut-shaped 
fullerene (Ceo) polymers exerts a strong influence on their low- frequency phonon states and their 
interactions with mobile electrons. The continuum approximation is employed to show the zone- 
folding of phonon dispersion curves, which leads to fast relaxation of a radial breathing mode in 
the ID Ceo polymers. We also formulate the electron-phonon interaction along the deformation po- 
tential theory, demonstrating that only a few set of electron and phonon modes yields a significant 
magnitude of the interaction relevant to the low-temperature physics of the system. The latter find- 
ing gives an important implication for the possible Peierls instability of the Ceo polymers suggested 
in the earlier experiment. 

PACS numbers: 72.80. Rj, 63.22.-m, 63.20.kd 
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I. INTRODUCTION 

Successful syntheses of fullerene (Cgo) molecules and 
related nanomaterials (T], Q have inspired tremen- 
dous interest in the community of nanotechnology and 
nanoscience. Experimental studies have shown that 
electron-beam irradiation to a Ceo film induces polymer- 
izations and produces one-dimensional (ID) Cgo poly- 
mers 01 • The polymers have peanut-like geometry, 
i.e., the tube structure with periodically modulated ra- 
dius along the axis ■ Such the peanut structure was 
confirmed by comparison between the infrared spectral 
measurement fsj and the first-principle density-functional 
calculations of the vibrational spectra '{I,, 9J. Recently, 
it was predicted that the periodic radius modulation 
strongly affects tr ansp ort properties of the carrier in the 
ID Cgo polymers |10[. Similar periodic corrugation of a 
two-dimensional conducting film was also shown to in- 
crease the electrical resistivity [llj. These findings are 
attributed to an effective electrostatic field caused by the 
periodic corrugation, which gives rise to a Fermi velocity 
shift and electron-electron scattering enhancement. 

The ID Cgo polymers exhibit semiconducting or metal- 
lic nature depending on the irradiation time interval 
This suggests that well-controlled irradiation time may 
lead the fabrication of quantum wires with desired con- 
ductivity. Of particular interest is the metallic Cgo poly- 
mers; in the metallic cases, nonequiliblium photo-excited 
carriers exhibit anomalously slow relaxation below 60K, 
indicating pseudogap formation in the electron energy 
band near the Fermi level [13] ■ Phenomenologically, this 
pseudogap is likely to result from the Peierls instability 
associated with a periodic distortion of the underlying 
carbon structure. However, its theoretical understanding 
has not been attained; in fact, it has remained unclear 
how the periodic radius modulation in the Ceo polymers 
affects the phononic excitations and electron-phonon cou- 
plings in the systems. 

We point out that the periodic geometry of the ID 



Ceo polymers will induce a pronounced shift in the low- 
frequency phonon spectrum where the phonon wave- 
length is comparable to the radius modulation period, 
i.e., larger enough than the carbon bond length. In the 
low-frequency range, continuum elastic approximation is 
a powerful means to analyze the phononic excitations. 
Furthermore, the approximation allows to use the defor- 
mation potential theory [13] that describes the electron- 
phonon interaction in a conducting medium. Combin- 
ing the two methods, therefore, we can formulate the ra- 
dius modulation effect on the phonon spectrum and the 
electron-phonon interaction in the ID Ceo polymers; the 
results obtained give a fundamental clue to explore novel 
properties of collective excitations in the system. 

In this paper, we study the low-frequency phonons and 
their coupling with electrons in the ID Ceo polymers by 
using a continuum model. A generalized eigenvalue equa- 
tion is formulated to obtain phonon dispersion relations 
that show zone-folding due to the periodic geometry. It is 
demonstrated that only a particular set of phonon modes, 
among many choices, yields marked enhancement in the 
electron-phonon coupling strength; these specific phonon 
modes are thought to play a prominent role in the pos- 
sible Peierls instability in the ID Ceo polymers. We also 
show that the zone-folding in the phonon spectrum en- 
courages the relaxation of radial breathing modes in the 
ID Ceo polymers, which is a consequence of the periodic 
radius modulation. 



In Sec. Ill A[ we derive an eigenvalue equation that is 
useful for calculating eigenmodes of general isotropic elas- 
tic medium with cylindrical symmetry. In Sec. Ill Bl we 
introduce a continuum elastic model of the ID Ceo poly- 
mers and suitable basis functions for an efficient eigen- 
mode analysis: Appropriate values of the elastic con- 
stants and geometric parameters of the model are pre- 
sented in Sec. Ill Cl Section Hill shows the zone- folding of 
phonon dispersion curves caused by the periodic geom- 
etry. Section. IIVI is devoted to formulate the electron- 
phonon interaction Hamiltonian on the basis of the de- 
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with the strain tensors [2l|: 
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FIG. 1. (Color online) Continuum model of a ID peanut- 
shaped Ceo polymer. The dashed square (red in color) indi- 
cates the unit cell. 



formation potential theory. Numerical results of the cou- 
pling strength are shown in Sec. |Vl and associated phys- 
ical implications are discussed in Sec. IVIl The paper is 
closed by conclusion given in Sec. IVIII 
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where di means the partial derivative with respect to 
the variables i = r,6,z. The integral in Eq. (jlj should 
be taken over the whole elastic medium. The first term 
of the integrand in Eq. ([1]) represents the kinetic energy 
density, and the other two terms account for the potential 
energy density with opposite sign. 

To the aim, we impose stress-free boundary conditions 
on the whole surface of the medium and follow the pro- 
cedure described in Ref . 22] . Expand Ui by a set of basis 
functions (/)Q,(r): 



(3) 



II. 



CONTINUUM ELASTIC MODEL 



In this section, we formulate a continuum elastic model 
of the ID Ceo polymers. Thus far, several contin- 
uum models have been developed to study low-frequency 
phonons in nanostructures such as GaAs/AlAs based 
nanowires ,15.1, pristine Cgo molecules [16| , and car- 
bon nanotubes |17H20j . Our model discussed below takes 
into account the effects of the periodic structure on the 
vibrational properties of the system, and the resultant 
vibrational eigenmodes will be applied for defining the 
electron-phonon interaction Hamiltonian in Sec. IIVI 



A. EIGENVALUE EQUATION 



where a labels a basis function. Preferred forms of (pai^) 
depend on the shape of the medium considered, and those 
for the present study will be given in Eq. Substi- 
tuting Eq. ([3]) into Eq. ([Ij and employing the variation 
method with respect to Aai, we obtain the generalized 
eigenvalue equation: 



(4) 



The matrix elements of 
ten as 



have a simple form writ- 



,dV. 



(5) 



Expressions of r(^i) (^jj) are somewhat complicated, thus 
summarized in Appendix \K\ 



The outline below is our derivation of a generalized 
eigenvalue equation [i.e., Eq. (U)] that allows to obtain 
the vibrational eigenmodes of an elastic medium. For 
general cylindrical objects, elastic vibration is governed 
by the Lagrangian: 



L = 



puj 



E 



Ui Ui 



2^ 



mEE^ 



dV. 



(1) 

Here, p is the mass density, w is the vibrational frequency, 
A and jj, are Lame coefficients. The complex-valued func- 
tion Ui denotes the ith component of the displacement 
vector u{r) at a point r, and u* is its complex conju- 
gate. Using the cylindrical coordinates, we can write 



B. HOLLOW TUBE WITH MODULATED 
RADIUS 

Next we restrict our attention to elastic thin hollow 
tubes with thickness 2h and length i, which we consider 
as a continuum model of the ID Cgo polymers. The tube 
axis points to the z direction, and the tube radius is 
spatially modulated as 



a a ( luz 

R[z) = i?o - - + - cos 

2 2 \ s 



(6) 



where Rq and Rq — a are the maximum and minimum 
radii, respectively. Since R{z) is periodic with a period 
s, we can define a unit cell by a portion of the tube within 
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the region z € [0, s] as depicted in Fig. [TJ For later use, 
we introduce a new variable ^(z) defined by 

i{z) = f w{z')dz', w{z) = ^Jl + (d,R)\ (7) 



Here, £,{z) measures the length of a geodesic curve [2J| 
between two points along the curved membrane repre- 
sented in FigUJ In terms of £,{z), a unit cell of the tube 
is a portion within the region ^ e [0, A] where A = ^(s). 

Once we specify the system's geometry as above, the 
volume integral in Eq. ^ is written by 



dV ^ I dz rdr / dO. 

V Jo JR(z)-h Jo 



(8) 



To eliminate the 
Ui{r,9,z) = Ui(r,t 



finite length effect, 
,z + L) with L = Ns, 



we suppose 
where N is 

the number of unit cells and = 800 in our calcula- 
tion. Under the condition, the wavenumber q along the 
z direction serves as an additional index (complemen- 
tary to the label a) that specifies phonon modes; here 
q = 2iTp/{NK) with -N/2 < p < N/2 {p is an integer). 
The basis function with q and a is defined by 



(9) 



where / = 0,1,2, 



m = 0,±1,±2,---, and G = 
27m/A {n = 0, ±1, ±2, • • • ). We used in Eq. (0) the ex- 
pression of e"?^ rather than e*''^, since phonons propa- 
gate along the ^ axis, not the z axis. Expression ^ im- 
plies that the label a represents the set of three integers: 
(l, m, n). 

It is noteworthy that (f>q^a with m and (fi* ^ with 
m) are orthogonal in the sense of the volume integral ([8|) . 
Therefore, individual phonon dispersion curves in the q- 
Lo plane are obtained by solving the eigenvalue equation 
(|4]) with m being fixed. 



C. NUMERICAL PARAMETERS 

In the phonon mode analysis, it is essential to deter- 
mine appropriate values of the Lame coefficients A, p, and 
the mass density p introduced in Eq.([T|). The Lame co- 
efficients are expressed by [2l| 



A 



{l + v){l-2v)' 



2{l + i^)' 



(10) 



where E and v are the Young modulus and the Pois- 
son ratio of the system, respectively. To our knowledge, 
there has been no attempt of determining E, v, and p 
of ID peanut-shaped Ceo polymers. Meanwhile, pris- 
tine Ceo molecules were suggested to have E = 1.06TPa, 
V — 0.145, and p — 2.27g/cm'^, whose values are simi- 
lar to those of graphene sheets 1^. We postulate that 
the three material constants of the ID Ceo polymer show 
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FIG. 2. (Color online) Phonon dispersion relations of hollow 
cylinders with periodically varying radius (a) and those with 
constant radius (b). The angular momentum index m is set 
to be m = for both plots. Six capitals (A-F) are assigned 
to the eigenmodes at gA = vr for each dispersion curve in (a) . 



the similar values to those noted above, motivated by the 
fact that adjacent Ceo molecules in the polymer are con- 
nected via the a bonds as same as adjacent carbon atoms 
within a pristine Ceo molecule. 

Given below is the list of geometric parameters we have 
used: Rq = 3.5A, s = 7.0A, and a = 1.5A, all of which 
were referred to the experimental characterization of the 
ID Ceo polymers discussed in Ref. fo*]. The tube thickness 
is fixed to be 2/i = 3.2 A as derived from the relation 
p — Mt/v, where Mt = 1.2 x lO^^'^kg is the total mass of 
60 carbon atoms and v ^ Att [{Rq + h)^ - {Rq - hf] /3 
is the volume of a hollow spherical shell with thickness 
2h that corresponds to the pristine Ceo atomic structure 
when it is smeared out. 

Before addressing the polymer's dispersion, we have 
checked the validity of our formulation by applying it 
to a straight hollow cylinder with no radius modulation. 
The exact solution of the vibrational eigenmode in the 
straight cylinder is g^("^s+<i^) multiplied by r-dependent 
Bessel functions (and their derivatives) [l^. In our 
formulation, the r-dependent Bessel terms are approx- 
imated by a set of polynomial functions, as seen from 
Eqs. ([3]) and ([9]). Despite this approximation, the com- 
puted phonon frequencies were in agreement with the 
exact solution to more than 99% accuracy when we set 



6. The same value of Imax has been used in the 



calculations of the modulated cylinders (i.e., the ID Ceo 
polymers), wherein 35 basis functions were employed to 
solve the eigenvalue equation in fact, we set < Z < 6 
and —2 < n < 2 for each fixed m. It should be noted 
that for very high frequency range (a; > 700 cm~^ in our 
study; see Sec. IIII[) . /max must be larger than our choice, 
but such the high frequency range is beyond the scope of 
the present analysis. 
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FIG. 3. (Color online) Dispersion relations for |m| — 1,2,3 

of hollow cylinders with: (a) periodically varying radius and mode E mode F 

(b) constant radius. 



III. DESPERSION RELATIONS AND MODE 
ANALYSIS 



Figure [DJa) shows the lowest six phonon dispersion 
curves for m = in the modulated tube. For compari- 
son, we plot in Fig. HJb) the lowest three curves for the 
straight tube with radius 3.5A; the six curves in Fig.[2lja) 
reduce to the three ones in Fig. [2Jb) if we set a — OA see 
Eq. ([6])]. In the low frequency range (w < 200 cm^^), 
the polymer's dispersion relations are almost identical 
to those of the straight tube; two almost-linear disper- 
sions are present. Among the two, the bottom curve 
belongs to a twisting mode (w = Ctq, Ct = \/TiJp — 
14.3km/s), and the upper to a longitudinal mode (w = 
ciq, ci — ^(A + 2/i)/p ~ 22.1km/s). Similar two linear 
branches were found in carbon nanotubes through con- 
tinuum [3, [T^ and atomistic calculations [24, 25], while 
the magnitude of frequencies deviates from the present 
results. An important observation in Fig.[2ja) andlSfb) is 
a significant difference in the dispersion profiles at inter- 
mediate frequencies (200cm~^ < w < 700cm~^) between 
the polymer and the straight tube. The periodic radius 
modulation causes zone-folding at qA = tt, resulting in 
additional branches in this frequency range. This im- 
plies that the radius modulation effect can be traced by 
investigating the phononic properties in the branches. 

Figures |31[a) and El^b) show the phonon dispersion re- 
lations of the modulated and straight tubes, respectively, 
for |m| = 1,2,3. Similarly to the case of to = 0, zone- 
folding at qA = TT and marked difference in the disper- 
sion profiles at the intermediate frequencies are observed. 
In the following, we specify all the vibrational states in 
the plots by three indexes: q, to, and the band index 
M — 1,2,---, which labels the A/th lowest curve for 
each TO. 

Figure [5] shows the displacement distribution in the 



FIG. 4. (Color online) Displacement distribution of the lowest 
six modes at qA = tt within the region z € [s/2, 5s/2]. The 
magnitude of displacement vectors is signified by the lengths 
and colors of arrows (red arrows are longer than blue ones in 
color). 



vibrational modes at gA = tt for each six curves presented 
in Fig.[2la). The profiles within the region z £ [s/2, 5s/2] 
are picked up for all the six modes that we labeled by 
A - F in order of increasing the eigenfrequency. The 
modes A and D correspond to twisting modes in which 
all vectors are parallel to the 9 direction. The other four 
are hybridized modes of longitudinal and radial breathing 
modes, where vectors point to directions in-between the 
z and r axes. 



IV. ELECTRON-PHONON INTERACTION 

In view of electron-phonon couplings, the hybridized 
modes of C and F depicted in Fig. 2] give an important 
consequence. To explain it, we formulate in this section 
the interaction Hamiltonian between electrons and low- 
frequency phonons along the deformation potential the- 
ory [ij]. The theory states that low- frequency phonons 
perturb electronic states through the occurrence of a de- 
formation potential represented by DV • u. Here, D is 
the deformation potential constant and u denotes the 
normalized displacement vector. The interaction Hamil- 
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tonian is given by 



= / vl/t(0, z) [D^ . u] -^[e, z)dS, (11) 



^c — ph 



where represents the annihilation (creation) oper- 

ator for electrons. The value of D lies within the range 
10-30 cV for various nanocarbon materials [IB, [13]) al- 
though its exact value for the ID Cgo polymer has not 
yet been examined; we thus set D = 20 eV in the present 
work. The surface integral in Eq. ()11|) is written as 



dS 



l'2Tr i-Ns 

de , 

Jo Jo 



dzR{z)w{z) 



(12) 



We now derive an alternative form of He-ph [i-e., 
Eq. (dH])] by expanding the field operators u and ^ in- 
troduced in Eq. (|lip . First, the zth component of the 
quantized phonon field u is expressed by 



J 



[bjuj.i + h\u*j^^ 



(13) 



where 5j(6j) is the phonon annihilation (creation) op- 
erator and J denotes the mode index of phonons, i.e., 
LOj ~ ijJm,q,M- The uj,i in Eq. (fT5|) is written as 



Uj,i = cj V^^-Va(r). 



(14) 
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FIG. 5. (Color online) Electron band structure and phonon 
emission process for several bands. The intra-band electron 
scatterings of (i), (ii), and (iii) within the electronic band 
labeled by (I7I, H) — {0,1), (0,2), and (1,1) are depicted, re- 
spectively. 



Finally, we arrive at the interaction Hamiltonian in the 
reciprocal-space representation by substituting Eqs. p3)) 
and ((ni) into Eq. ([TT|) : 



Following the quantization procedure [28|, cj should be 
determined by carrying out the integral: 



J Uj, ■ ujdV = 



2pujj 



'J.J' 



(15) 



Consequently, we obtain the normalization constants: 



1 



cj 



(16) 



where A'^'^^ is the eigenvector corresponding to uij, E is 
the matrix defined in Eq. (U), is the volume of the 
system. 

Next, the field operator for the electron is written as 



(17) 



where is the annihilation operator of the electron 

with spin a, 7 is the angular quantum number, k is the 
wavenumber, and H is the band index. Explicit forms 
of the electron's eigenfunctions ipl k ^) derived in 



Appendix |B] The electronic band structure eiy as a 
function of k is obtained by numerical diagonalizations of 
the Schrodinger equation (|B1I) . The lowest few electronic 

dispersion curves e-y (fc) for several 7 and H are drawn in 
Fig. [S] Zone-folding at kA = ±7r is due to the curvature- 
induced electrostatic potential with the same periodicity 
of the radius modulation [l^, . 



where the labels J — (to, q, M) and Pj = (7^, kj, Hj) in- 
dicate the phonon and electron modes, respectively. The 
p. 

quantity g j" ' describes the strength of the electron- 
phonon interaction through which electrons are scattered 
from initial states [j — i) to final ones (j = /); several 
scattering processes are illustrated in Fig. [5] We can 
prove that 



Rq ' '-I'G.O 



(19) 



where 



K > 



Pi ,Pf .771, q 



l,G,p 



E^- 



7i+m,7/0fei+q-fe/,27ri/A 



G',G" 



X / dzw 
Jo 



1+p 



Rn 



J(2ttI/ K-G+G' ~G")(,(z) 



(20) 



and C^^L g is a Fourier expansion coefficient of ip j, as 

given by Eq. (jB3p . In Eq. (PH)) , the integer / is introduced 
to let ki and fc/ be in the first Brillouin zone, i.e., kjK e 
[-7r,7r]. 
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FIG. 6. (Color online) Coupling strength \g\^ in the case of FIG. 7. (Color online) (a) The fci-dependence of for j7| = 

7 = 0. The dependences on the electron's initial wavenumber 1. (b) Rapid decrease in jgp of M = 6 at fc^A — -k 12 with the 

ki are displayed, while the phonon mode parameters m = electron band index H. 
and gA = TT are fixed. 



V. COUPLING STRENGTH 

The coupling strength for various conditions 

of m,q,M) and Pj{= jj,kj,Hj) with j = i or f 
has been systematically calculated, only a part of which 
are presented in Figs. [6] and [T] In all the figures, we set 
m = 0, gA = TT, 7i = 7/(= 7), and Hi — Hf{= H), since 
the intra-band electron scatterings that satisfy these con- 
ditions (see Fig. [S|) were found to yield significantly large 
values of \g/' ^p. Hereafter, the suffixes m g j" ' are 
omitted for simple notation. 

First, we look into jgp for 7 = and iJ = 1, in 
which only the process (i) depicted in Fig. [5] is rele- 
vant. Figure ini[a) shows jgp as a function of the initial 
wavenumber ki for four different Ms, i.e.. Pi = (0, hi, 1), 
Pj = (0, fcj + tt/A, 1), and J = (0,7r/A,M). The con- 
tributions from M = 1,4 to \g\'^ vanish at every ki, be- 
cause these twisting modes give V • m = 0. For almost 
all ki, with M = 3 takes the largest value, while 
|gp with M = 6 competes within limited ranges around 
fc^A ~ ±7r/2. The primary contributions to \g\'^ from the 
phonon modes of M = 3, 6 can be explained by their lon- 
gitudinal nature. As shown in Fig. |4l the displacement 
vectors in the two modes (C and F) are mostly tangential 
to the modulated cylindrical surface, thus yielding large 
values of V • M. This results in the enhancement of the 
electron-phonon coupling as follows from the expression 
of Tie-ph given in Eq. ([T^. 

The profile of \g\'^ drastically changes when electrons 
stay at upper bands of > 2. Figure ^h) shows the 
contribution to |gp from the process (ii) in Fig. [SJ i.e., 
H — 2 and 7 = 0. In this case, I5P is strongly depressed 
except for a few peaks at fc^A = 0, ±7r. If we increase H 
further, the peaks diminish rapidly and jgp disappears 



for all ki. This disappearance is in contrast with |gp 
for H ^ 1, where every ki give finite values of \g\'^. We 
thus conclude that, in the modulated cylinders, signifi- 
cant values of \g\^ can be obtained only when the coupled 
electrons belong to the lowest band of = 1. 

Next, we pay attention to the case of 7 7^ 0. Fig- 
ure [Tl^a) shows |(?ps for I7I = 1 and H = 1 [the process 
(iii) in Fig. [5], in which |gp with M = 6 is the largest for 
all ki and that with M = 3 shows a modest value. We 
have confirmed that as H increases, every curves of jgp 
go downward, which result in profiles of |gp analogous 
to Fig. [6jb). Similar depression of \g\'^ with increasing 
H have been observed for I7I > 2, as shown in Fig. [Tl^b). 
We can say that for arbitrary 7, \g\'^ for A/ = 3 or 6 has 
significant value only when both two coupled electrons 
reside at a small H branch. 



VI. PHYSICAL CONSEQUENCES 

A. DOMINANT MODE COUPLING 

Our findings of \g\'^ in Sec IVl give an important sugges- 
tion for addressing the property of the possible Peierls 
transition in the ID Cgo polymers. As is well known, the 
Peierls transition temperature Tc in general ID systems is 
dependent on the coupling strength |oP; the larger value 
|(7p arises, the higher Tc we observe [30|. It follows from 
Sec. |V] that \g\'^ becomes significant only when the elec- 
trons lie at a small H branch and a specific phonon mode 
of Af = 3 or 6 is involved. All the rest combination of 
electronic and phononic modes give slight contribution to 
the Peierls transition in the ID Cgo polymers. In this con- 
text, we have succeeded to determine what specific modes 
(electron and phonon) in the ID Cgo polymers should be 
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FIG. 8. (Color online) Allowed relaxation process of RBMs 
for (a) the ID Ceo polymer and (b) a straight tube. The L-, T- 
, and F-processes suggest the RBM's relaxation in which the 
post-decay modes are the longitudinal, twisting, and flexural 
modes, respectively. 



responsible for the Peierls transition endowed with a fea- 
sible value of Tc. Quantitative evaluation of Tc for the ID 
Cgo polymers is possible via the standard Green's func- 
tion technique [31j incorporated with our formulation of 
|gp given in Eqs. ([T9)) and pO|. The calculated results 
of Tc as well as the comparison with experiments will be 
discussed elsewhere. 

We emphasize that the present formulation of jgp is 
broadly applicable to the study of the radius modula- 
tion effects on the physical properties of ID Ceo poly- 
mers. The application should include the electrical resis- 
tivity analyses. In general, the low-temperature resistiv- 
ity of nanomaterials is dependent on temperature and the 
Fermi energy, and the dependences are governed largely 
by the magnitude of jgp. Hence, our formulation of |gp 
will be helpful to estimate the low-temperature resistiv- 
ity of the ID Ceo polymers that remains untouched to 
date. In addition, artificial control of the Fermi level by 
chemical doping or bias voltage application may enable 
to obtain a desired magnitude of leading to an efB- 
cient means of developing next-generation nanodevices. 

Besides the modulated cylinder problem, an extension 
of the present theory to three dimensional counterparts 
(e.g., the Mackay crystal [s^l and the zeolite-templated 



carbon [33|) is quite interesting. The Mackay crystal is 
a three-dimensional carbon network; three infinite tubes 
intersect smoothly in the unit cell, resulting in a triply 
periodic minimal surface. Due to the intriguing geome- 
try and functionality [s^l, it has attracted interest since 
the first hypothesis. We speculate that in the Mackay 
crystal, only specific phonon modes cause large deforma- 
tion potential similarly to the present system, since each 
constituent tube in the three dimensional minimal sur- 
face can be regarded as a radius-modulated tube along 
the axial direction. Electronic structures of the periodic 
minimal surface have been already clarified [3514371 , and 
thus employing the results to obtain the electron-phonon 
interaction Hamiltonian will yield a conclusion of this is- 
sue. 



B. RADIAL BREATHING MODE RELAXATION 

Another possible consequence of the periodic radius 
modulation in ID Ceo polymers is a considerable decrease 
in the relaxation time of radial breathing modes (RBMs) 
compared with that in carbon nanotubes [38|. We have 
numerically found that the phononic mode labeled by 
(m, q, M) = (0, 0, 3) is a RBM, which is active on the reso- 
nant Raman scattering. The lifetime of a RBM is mainly 
determined by phonon-phonon scatterings which lead to 
a phonon decay |38)]. When the number of the allowed 
decay process is large, the corresponding phonon mode 
has a short life. This indeed holds true for the RBM of 
the ID Ceo polymers, as demonstrated below. 

Figures HJa) and Mjo) illustrate the RBM decay pro- 
cesses in the ID Ceo polymer and the straight tube, re- 
spectively. Allowed decay processes are represented by 
slanted arrows that extend from the RBM (indicated 
by u!fi) to low- lying acoustic modes having the eigen- 
frequency = ujr/2: The post-decay states locate at the 
intersections of the horizontal dotted line (= ujr/2) and 
dispersion curves. Note that the decay process is con- 
strained by the conservation laws of energy, wavenum- 
ber, and angular momentum. Therefore, a RBM hav- 
ing the eigenfrequency ujji decays into a pair of acous- 
tic phonons with ujm.q.M and uj-m.-q+G.M that satisfy 
Wfl = uJm,q,M + ^-ni-q+GM; here, G = Q for normal 
processes and G ^ for umklapp processes. It is found 
that in the straight tube (w^ = 350cm~^), there are 
four decay processes: L-, T-, and doubly-degenerate F- 
processes (i.e., wi^^^ + and 4- W-i,g,i), 

where L, T, and F mean that the post-decay modes are 
the longitudinal, twisting, and flexural modes, respec- 
tively. All the processes belong to the normal process 
(i.e., G = 0). On the other hand, the RBM of the poly- 
mer {loji — 440cm~^) can decay through both the nor- 
mal and umklapp processes (G ^ 0) involving L- and 
T-processes. Due to many choices of G available, the 
number of decay process in the polymer is much larger 
than that in the straight tube, although the F-process in 
the polymer is forbidden as the dispersion curves with 
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|m| = 1 has the frequency gap around ujr/2. In other 
words, the zone-folding in the phonon dispersion results 
in a considerable increase in the number of the decay 
process in the polymer case, which leads to a broadening 
of the RBM peak linewidth. It is known that at room 
temperature, the linewidth of the RBM peak in carbon 
nanotubes is ~ lcm~^ independently of the tube radius 
[ssj . Hence, that in the ID Cgo polymer is expected to 
be larger than lcm~^ as a result of the periodic radius 
modulation. The theoretical estimation of the linewidth 
will be carried out in our future work. 



VII. CONCLUSION 

We have investigated the low-frequency phonon states 
and their coupling with electrons in the ID peanut- 
shaped Cgo polymers along the continuum approxima- 
tion. It has been found that the phonon dispersion rela- 
tions exhibit the zone-folding effect due to the periodic 
radial modulation. The electron-phonon interaction was 
formulated in line with the deformation potential theory, 
and the obtained formula for the coupling strength allows 
us to identify which phonon modes should be responsible 
for the Peierls transition observed in the prior experi- 
ment. The present formulation for the coupling strength 

I 

For r(^j) with i = r, we have 

^(I3r),{az) 

For T{j3i)^{ai) with i 

^(130): 

For r(^i)_(„i) with i 

^{Pz 
Y{l3z 

Here, we have used the following expression 

$p,± =9,,^p±^. [p = a,P] (A4) 



has a wide applicability for estimating physical quantities 
of the systems, thus playing a crucial role in predicting 
untouched behaviors governed by the electron-phonon in- 
teraction. 
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Appendix A: EXPRESSIONS OF r(^i)(^;) 

The matrix elements r(^,i) (^j;) introduced in Eq. (|4]) 
are represented as follows. 



I 

The volume integral and the basis function ipadi) are de- 
fined in Eqs. ([U and ([9]) in Sec. IIIBl respectively. 



'pV 



J ^p,+d,KdV + d,cj)}dr<l>c.dV. 



(Al) 



A r 1 



(ar) ~ 

(az) = 
= -2, 
),(ar) ■ 

)Ac^z) ■ 



pV J r 
A /■ 1 



-de(b%^o..+dV 



de(j)*fjde(j)adV 



pV 



pV 



•a -r ^13 



-^l,_-decj)^ ) dV, 



^dg(j)ld0<j)c. + + dz<l)pdz(t)a dV, 



A ^1 

pV J r 



p fl 



— I -Mpd.'l^cdV + J -d,4>}deKdV. 



^ j d,c^l^^^+dV + j {dA}dAc)dV, 
j ^d^4>*pde(j)adV + ^ j ^de(t>}dz4>adV, 



dr4>*pdr4>a + ^^e^^^e^a + 2d,_(t)*pdz<j)a 1 dV. 



(A2) 



(A3) 
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Appendix B: ELECTRON EIGENSTATES 



is written as [l^, l29| 



2m* 



(Bl) 

Here, m* is the effective mass and U^{z) is an effective 
potential expressed by 



In this Appendix, we derive electronic states in the 
cylindrical tube with varying radius. We have assumed 
that the tube thickness is small enough to suppress elec- 
tron's excitation in the normal direction. The same 
assumption works well for treating the interaction be- 
tween electrons in both ID Cgo polymers and two- 



dimensional curved thin films Due to the rota- 



U,{z) 



i?2 



I? - RdlR 
2Rw^ 



(B2) 



Since R{z) and w{z) are periodic as shown by Eqs. ^ 
and d?]), U^{z) is also periodic, i.e., U-y{z) = U-y{z + s). 
Thus, it follows from Bloch's theorem that 

„i^0 



/27r 



1 



^/TTar 



tional symmetry of the tube, the wavefunction becomes 
^p{9,z) = {e^^^ / ^/2tt) Xi{^) with the electronic angular 
quantum number 7. The Schrodinger equation for x-i{^) where Ge = 27rne/A [n,, — 0, ±1, ±2, ■ 



(B3) 
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